X-ray amplification from a Raman Free Electron Laser 
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We demonstrate that a mm-scale free electron laser can operate in the X-ray range, in the in- 
teraction between a moderately relativistic electron bunch, and a transverse high intensity optical 
lattice. The corrugated light-induced ponderomotive potential acts simultaneously as a guide and as 
a low-frequency wiggler, triggering stimulated Raman scattering. The gain law in the small signal 
regime is derived in a fluid approach, and confirmed from Particle-In-Cell simulations. We describe 
the nature of bunching, and discuss the saturation properties. The resulting all-optical Raman 
X-ray laser opens perspectives for ultra-compact coherent light sources up to the hard X-ray range. 
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The advent of fully coherent light sources in the X- 
ray range promises to be the next revolution in X- 
ray science, leading to many scientific, industrial and 
health applications. Large scale X-ray free electron laser 
(XFEL) projects have been launched, and start supply- 
ing high brightness beams for novel physics experiments 
[1]. However, these large scale infrastructure cannot al- 
low for widespread dissemination of the XFEL technolo- 
gies, which would require to shorten very significantly the 
length of the linear accelerators and magnetic undulators. 
This size constraints also prevent to reach the hard X-ray 
range, which would open entirely new perspectives. 

Few alternative strategies have been proposed to sup- 
ply compact X-ray free electron lasers, based on substitu- 
tion of the LIN AC by laser- wakefield acceleration [2, 3], 
use of original compact undulators as an ion channel [4] , 
or a counter-propagating laser [5-7]. Laser undulators 
offer indeed a key advantage : laser wavelengths at the 
/am level allow one to reach X-ray photon energies with 
moderately relativistic electrons, of kinetic energies of 
few tens of MeV only. In all cases, the XFEL is expected 
to operate in the conventional regime of stimulated in- 
verse Compton scattering, imposing severe limitations on 
electron emittance, transport, kinetic energy spread, and 
laser uniformity for laser undulators [8], thus hindering 
prospects of experimental demonstrations. 




FIG. 1. X-ray Raman scattering geometry in a reference sys- 
tem moving with the electron bunch. 

A conceptually different new scheme considers a rel- 



ativistic electron bunch injected into the overlap region 
between two transversally incident, counter-propagating 
intense lasers beams [9]. The setup is depicted in Fig. 1, 
directly in the reference frame of the electron bunch. The 
interference between the laser beams forms an optical 
lattice, and induces a spatially corrugated ponderomo- 
tive potential for the incident electrons that is trapping 
them transversely. The electron dynamics then consists 
of high frequency oscillations induced by the two lasers 
along the laser polarization direction, and of low fre- 
quency oscillations along the interference direction, simi- 
lar to betatron oscillations, with a characteristic angular 
frequency Q. Light is hence scattered at the betatron fre- 
quency, and on the Stokes and ant i- Stokes lines around 
the laser frequency. In the laboratory frame, this scat- 
tering is Doppler up-shifted by 27 2 , where 7 is the elec- 
tron Lorentz factor. Scattering is spontaneous as long 
as the electron motions are uncorrelated; however, elec- 
trons may also exhibit a collective low-frequency oscilla- 
tory behaviour, so that we can expect a stimulated Ra- 
man instability and coherent emission of X-ray radiation 
in the forward direction. This Raman-type scattering 
should be distinguished from the known Raman instabili- 
ties in conventional long wavelength Free Electron Lasers, 
where the system oscillations are the Langmuir plasma 
waves [10]. This new Raman instability dominates if the 
bounce frequency Q is greater than the electron beam 
plasma frequency u p . 

The aim of the present Letter is to investigate the X- 
ray amplification process in this trapped mode on the ba- 
sis of particle-in-cell numerical simulations, and to derive 
its main scaling laws on small-signal gain and saturation, 
from an analytical hydrodynamic model. 

Let us consider a relativistic electron bunch moving 
along a z-axis, with an average Lorentz factor up to a 
hundred typically, and number density n® : incident onto 
the optical lattice resulting from overlapping twin laser 
beams, interfering along the x-axis. The light lattice in- 
tensity is considered as sub-relativistic, typically in the 
range 10 16 to 10 18 W/cm 2 for near infrared lasers. 

The process is conveniently described in the frame 
moving with the electron beam velocity = c/3& < 1, 
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for which the incident laser and scattered field frequen- 
cies are similar, and electron motion can be treated 
non-relativistically. The twin lasers appear in that 
frame at oblique incidence, with an angle defined as 
fco_i_/fco|| = VtA> ^ 1> see fig- 1- Sub-relativistic lat- 
tice intensity defines a moderate vector potential, which 
is a Lorentz invariant and in dimensionless units reads 
a = 0.85 x 10- 9 X l Vf, where I (W/cm 2 ) is the laser 
intensity and X l (/im) the laser wavelength. With the 
superscript we denote the parameters in the labora- 
tory reference frame. 

In that frame, electrons first enter a ramp-up region, 
followed by an optical lattice of a constant intensity; in 
the bunch frame, this process is described as a relatively 
slow switch-on of the external electromagnetic field. The 
normalized vector potential of the incident laser irradia- 
tion reads as 

&l = cli + &2 — 2ao(£) sin(/co±x) cos(co>o^ + k \\z) •> (1) 

resulting into a ponderomotive potential U p : 

U p (x, t) = m e c 2 ao(t) 2 sin 2 (ko±x) . 

Near the bottom of the potential, electrons oscillate with 
the bounce frequency Q(t) = V2ao(t)koj_c. 

Field amplitude grows during a time ti n j up to a$. 
Electrons with transverse velocities, v± < V^&oc, may 
be trapped in the ponderomotive potential and oscillate 
along the closed phase trajectories in (x,p x ) -plane. Par- 
ticles with lower initial velocities are trapped before U n j , 
and further growth of potential adiabatically compresses 
their phase trajectories, thus increasing the velocity am- 
plitude and electron density at the beam axis [11]. The 
maximum electron density after injection can be esti- 
mated as: 

n e = n° e (V2a /5f3 0± ). (2) 

Trapped electrons are driven at high frequency by the 
lasers, and oscillate transversely in the optical lattice po- 
tential, akin to the betatron oscillations in plasma ion 
channels [12, 13], with a harmonic law x e = x m cos(^t), 
where Vt = V^^o^o^c. Radiation is emitted in two fre- 
quency ranges : at the betatron frequency [14, 15]; and 
symmetrically around the laser frequency, at the Stokes 
ujo — Q and anti-Stokes ujo + £1 frequencies [9]. 

We have developed a numerical tool, dubbed EWOK, 
to investigate scattering process in the bunch frame. 
The two-dimensional code is based on a general particle- 
in-cell approach [16]; it factors out the high-frequency, 
laser-induced oscillations of electrons, solving rela- 
tivists motion equations for macro-particles in elec- 
trostatic and ponderomotive potentials. The scat- 
tered wave is described in the envelope approximation 
a s = aexp(— iuj^t + ik^z) as a slowly oscillating ampli- 
tude a. This allows us to reduce to first order the prop- 
agation equation along z, and to solve it over x by a 
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FIG. 2. Evolution of vector potential at the center of right 
boundary (z,x) = (Z z ,0). The exponential fit T/cjo = 1.8 • 
10 -3 is shown by a red dotted line. 

Fourier transform d 2 a^ —k 2 x a\ Zx . The simulation do- 
main is a rectangular box, whose dimensions l ZjX de- 
fine a discrete number of electromagnetic eigenmodes. 
Boundary conditions are periodic for both electromag- 
netic waves, while for the particles, boundaries are peri- 
odic along z-axis and absorbing along x-direction. Thus, 
the beat between the laser and the scattered wave must 
be periodic, which imposes constraints on the wave- 
vectors (k s — ho) of the scattered wave, and (k s + & |i) °f 
the longitudinal ponderomotive force. For a plane scat- 
tered wave at the Stokes frequency, the corresponding 
periodicity conditions read: 

fc fl = A*(l-WiA /Z*) (3) 
lz/*o = N 2 /(l + i3 b ). 

Ni and N 2 being integers. These conditions lead to ex- 
ceedingly long simulation domains for large 7 values. We 
perform therefore simulations with a relatively low 7, and 
exploit the analytical laws to rescale the results to exper- 
imentally relevant parameters. 

We have used EWOK to simulate amplification on 
the Stokes mode, satisfying the resonance condition 
ujo — uj s = 1), with physical parameters summarized in 
table I, and on the basis of 5.2 • 10 6 macro-particles. 

Figure 2 presents the time evolution of the scattered 
field amplitude on the right boundary. In the interval 
100 < t < 500Ao/c, one can see an exponential growth 
of the signal, followed by saturation at the amplitude 
(a s )/ao ~ 1.4 • 10 -3 . The exponential growth rate is 
estimated as T/ujo = 1.9 • 10 -3 , the fit being plotted as a 
dotted line. For 1 /im laser wavelength, this corresponds 
in the laboratory frame to a gain length in intensity of 
° nl y L gain = 84 /im. 

The laser electric field changes its sign across the zero 
lines of ponderomotive potential wells. The beat be- 
tween the laser and the scattered wave is hence in phase 
opposition on the two sides of a well, resulting in an 



3 




65 66 67 68 69 70 



FIG. 3. Snapshot of electron density distribution at satura- 
tion time 600Ao/c. 



TABLE I. Interaction parameters 



laboratory moving beam 



electron beam 


current 


8 kA 


density 


2 • 10" 5 n c 


width 


3/xm 


width 


3/im 


duration 


10 fs 


length 


30 /im 




io- 4 


SP\\o 


io- 4 


emittance 


0.3 mm-mrad 




0.1 


energy 


4.7 MeV 


7 


9.4 


optical lattice 


\ l 
A o 


1 /xm 


A 


106 nm 


intensity 


1.6- 10 16 W/cm 2 


ao 


0.11 


ramp 


200 fim 


Tinj 


200 X l /c 


scattered light 


Ai 


5.6 nm 


0J S 


0.985 u 



ant i- symmetric bunching structure. This is illustrated 
in fig. 3, showing the completely trapped electron distri- 
bution at saturation. A quasi-sinusoidal bunching shape 
can be noticed, in sharp contrast with the series of par- 
allel micro-bunches of conventional free electron lasers. 

These numerical findings can be completed by an an- 
alytical analysis; while several theoretical approaches 
are possible, we choose to put forward a hydrodynamic 
model, that yields relevant estimates for the scaling laws 
within a simple mathematical framework. We propose to 
model amplification analytically by describing the collec- 
tive electron behaviour n(x,z,i) as an eigenmode of the 
light potential well. In a simplified approach we neglect 
the electrostatic potential of electrons and the longitu- 
dinal temperature, and consider the transverse tempera- 
ture not to exceed the trapping limit. Electron motion 
follows hydrodynamic equations, coupled to the propa- 
gation equation for the light vector potential : 



d t u x + V x P xx /m e n e + V ' X U P = (4b) 
d t u z + V z U p = (4c) 

((d t - iuJo) 2 - c 2 V 2 ) a = -(47re 2 /m e )n e a (4d) 

where u is the electron fluid velocity. Considering the 
slowly oscillating vector potential a defined earlier, the 
ponderomotive potential reads U p = m e (c/2) 2 |a| 2 . Note 
that a in eq. (4) is the total field, including the scattered 
and laser fields. 

The beam is trapped along x, so it follows the adia- 
batic equation of state with the electron pressure reading 
Pxx — vnv\v? / 3no, where v± = c5f3_\_o is defined by the 
transverse velocity spread. Considering equilibrium be- 
tween lattice potential and electron pressure, we find the 
transverse density distribution of the background state 

n e riQ\/l - £ 2 , 

where ( = x/L is a coordinate normalized to the beam 
width L = v±/Q, and no = n e (0) is a maximum electron 
density. 




FIG. 4. Spatial-spectral distribution of electron density per- 
turbations. The predicted coupling shape g(£) is shown in 
white curve for L = Ao and the red dashed line is a slice of 
the spatial-spectral diagram along x-axis at uj = Q. 

In the linear approach, we consider the collective elec- 
tron modes as a first-order perturbation of the beam po- 
sition Ajc, defined as £ = (x + A x )/L. Without a signal 
wave there is no longitudinal force d z U p , and from eq. (4) 
it follows: 

which describes the collective beam oscillations with a 
frequency Q. 

The small-amplitude signal wave a s is taken as a per- 
turbation to the laser vector potential, inducing to first 
order a ponderomotive potential = m e (c/2) 2 a s a^ + 
c.c. Assuming a plane scattered wave - thus, neglecting 
diffraction -, we linearize eq. (4) and turn to the Fourier 
domain, to obtain the dispersion equation: 



d t n e + V • (n e u) = 



(4a) 



(uj 2 - tt 2 ) ((co - uu ) 2 - (k z - k 0ll ) 2 c 2 ) = au$ . (5) 
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This equation describes electron and electro- 
magnetic modes with a coupling coefficient 
a = G(uj p aok z cko_\_L/ujQ) 2 , where G is an overlap 
integral resulting from an inhomogeneous profile of the 
laser field oc x and the beam density ex y/l — x 2 /L 2 . 
Thus, the model predicts that amplification is localized 
mainly in the regions around the coupling maxima at 
£ = iy/2/3, with widths £ ~ 0.5. This analytical 
prediction was checked numerically. Figure 4 shows 
the spectra of electron density perturbations (abscissa) 
as a function of the transverse position for an 
arbitrary fixed position z. The spectra are centered at 
the betatron frequency Vt = 0.016 c^o- The transverse 
distribution consists of two symmetric maxima, as 
predicted by the hydrodynamic model. The dashed line 
is obtained from the simulation, while the white curve 
shows the transverse coupling function for a trapped 
beam half- width L = X . The numerical shape is slightly 
broader, but the general structure of the emitting zones 
is well reproduced. 

Solution to eq. (5) for the Stokes mode is a complex fre- 
quency with a real part uo ~ Q and a positive imaginary 
part T, defining the amplification growth rate. Assuming 
a weak coupling condition a <C (ft/uo) 3 , the growth rate 
is: 

T/uo =0.5(au /n) 1/2 . (6) 

Estimating G c± 0.2, and using parameters applied in 
numerical modelling (table I), we obtain a growth rate 
T/uo = 1.7- 10- 3 . 

The exponential growth of the signal saturates when 
the electron beam gets trapped in the longitudinal po- 
tential |a s a|J, resulting a full particle bunching n e ~ 
noexp(— iQt) + c.c. From eq. (4d), the maximal signal 
amplitude can be deduced as : 

(a s )/a ~ (k ±L)ujp/(Slujo) • (?) 

For the chosen parameters, the signal saturation occurs 
at (a s )/ao = 1.7 • 10 -3 . These analytic estimates of 
growth rate and maximum of the scattered field are in 
good agreement with results of numerical modelling. 

In the laboratory frame, the growth rate reads T l /uj 1 q — 
T/2o;o , and can be expressed as a function of initial 
parameters: 

r74 = 0.1 (Kh) (Ao/o-j.) y/ao J I J a (8) 

where K = ko±L is a numerical factor defined by the 
filling of potential channel by electrons and J a = 17 kA 
is the Alfven current. In normal trapping conditions, 
K ~ 1. 

These scaling laws provide estimates of the gain and 
saturation level from a Raman XFEL laser derived from 
laser-wakefield acceleration. Let us consider parameters 
close to ones achievable in the experiments [17] : an elec- 
tron bunch of 50 MeV, with a peak current of 14 kA, 



and 1 mm-mrad emittance. Twin undulator lasers are 
at a wavelength of 800 nm, with an intensity of 3 • 10 1T 
W/cm 2 , corresponding to a normalized vector potential 
of 0.4. The output X-ray photons have an energy of 30 
keV; the gain in intensity extrapolated from our numer- 
ical results is 15 cm -1 . Assuming a transverse size of 3 
/im, an X-ray duration of 10 fs, results in an X-ray energy 
at fiJ level, or 10 8 — 10 9 photons. 

Note that photon emission results in an electron recoil. 
In the electron rest frame, the recoil can be observed nu- 
merically as a collective backwards motion along z-axis, 
which starts near the saturation stage. Since emission is 
a quantum process, in the case of high photon energies 
the recoil has a discrete nature, and a quantum descrip- 
tion may become appropriate. A significance of quantum 
effects is defined by a quantum-recoil parameter [18], as 
a ratio of the recoil induced detuning and the gain band- 
width. It reads in our model q' x = ( / y e huo)/(m e c 2 T/ujo). 
Parameters in the present examples are restricted to 
q x < 1, where a classical description is valid. 

In summary, we have shown that a relativistic electron 
bunch, injected into and guided by a high intensity opti- 
cal lattice, triggers a forward Raman instability resulting 
in a rapid amplification of a coherent X-ray beam. If cou- 
pled to laser-wakefield acceleration, this process promises 
to yield hard X-ray laser beams. Raman processes of- 
ten exhibit an intrinsic robustness with respect to ran- 
dom variations of the driving parameters [19], such as 
the laser intensity, which may be a key to an experimen- 
tal demonstration. The Raman X-ray free electron laser 
should therefore be further studied as an opportunity to 
supply ultra-compact coherent X-ray sources, up to the 
hard X-ray range. 
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